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Charge symmetry breaking from a chiral extrapolation of moments of quark 

distribution functions 
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We present a determination, from lattice QCD, of charge symmetry violation in the spin- 
independent and spin-dependent parton distribution functions of the nucleon. This is done by 
chirally extrapolating recent QCDSF/UKQCD Collaboration lattice simulations of the first several 
Mellin moments of the parton distribution functions of octet baryons to the physical point. We 
find small chiral corrections for the polarized moments, while the corrections are quantitatively 
significant in the unpolarized case. 
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I. INTRODUCTION 



Charge symmetry refers to the equivalence of u quarks 
in the proton and d quarks in the neutron, and vice- versa. 
Precisely, it is the invariance of the strong interaction un- 
der a rotation of 180° about the '2'-axis in isospin space. 
At low energies, charge symmetry is obeyed to a preci- 
sion of order 1% [TJ. It would be natural to expect that 
partonic charge symmetry should hold to a similar ex- 
tent. Traditionally, this expectation has been applied to 
parton phenomenology [2,3, and the assumption of good 
charge symmetry has been used to reduce the number of 
independent quark distribution functions by a factor of 
two. 



Recently, charge symmetry violating (CSV) effects 
have been included in phcnomcnological parton distri- 
bution functions for the first time [3] , and theoretical es- 
timates of the size of such effects have been made [3 E] . 
Experimental upper limits on partonic CSV are in the 
range 5-10% CSV of this magnitude would produce 

important effects in tests of physics beyond the standard 
model, for example in neutrino deep inelastic scattering 
experiments [7]. 

The first two Mellin moments of the spin-dependent 
quark distribution functions of the octet baryons, and 
the second spin-independent Mellin moment, have re- 
cently been determined from Nf = 2 + 1 lattice simula- 
tions by the QCDSF/UKQCD Collaboration 8, 9 . The 
first analysis of this lattice data used a linear flavour ex- 
pansion about the simulation SU(3) symmetric point to 
extract values for the charge symmetry violating distri- 
butions 8 . Using chiral extrapolation formulae for the 
Mellin moments of quark distribution functions [TDlU6| . 
recently extended to include CSV effects [T7] , we improve 
on the original analysis by extrapolating the lattice re- 
sults to the physical point. We find that chiral physics 
generates small corrections to the parton CSV terms. 



II. METHOD 

In terms of quark distributions, charge symmetry im- 
plies 

u p (x,Q 2 ) = d n {x,Q 2 ), d?{x,Q 2 )=u n {x,Q 2 ), (1) 

with analogous relations for the antiquark distributions. 
A measure of the size of the violation of charge symmetry 
is given by the 'CSV parton distributions', defined in 
terms of the Mellin moments as 



Su m± = dxx m (u p± (x)-d n± (x)) 
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dxx m {d p± {x)-u n± {x)) 



o 



Of 



(2) 
(3) 

(4) 
(5) 



for the spin-independent distributions, with analogous 
expressions for the spin-dependent case. Here, the plus 
(minus) superscripts indicate C-even (C-odd) distribu- 
tions: 



Q ± {x) = q(x)±q(x). 



(6) 



The Mellin moments accessible to lattice simulations 
alternate between C-even and -odd moments with in- 
creasing m; the (m — l)th spin-independent (SI) and mth 
spin-dependent (SD) lattice moments are defined as 



„m-l\B 



(x m )i q 



dxx 
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{q a {x) + {-l) m r{x)) 



(7) 



dxx m {Aq B (x) + {-l) m Aq B {x)). (8) 



Recent lattice simulations from the QCDSF/UKQCD 
Collaboration jS] give results for the first two SD and first 
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SI lattice moments. As these Nf = 2 + 1 lattice simula- 
tions use degenerate light quarks, the CSV terms cannot 
be directly evaluated from the simulation results using 
Eqs. |2]) and ([3| (as this would give zero in each case). 
The problem can, however, be approached indirectly. 

The original analysis of the QCDSF/UKQCD Collab- 
oration lattice data used a linear flavour expansion at 
the simulation SU(3) symmetric point to estimate the 
CSV terms [5]. That is, the CSV terms were expressed 
in terms of hyperon moments as 



Su = ms ( — 



dm u 



dm d 



0{m 2 ) 



(9) 



where m$ = (rn d — m u ), and Sd may be written in a 
similar way. The equivalence of the u and d quarks in 
the simulations, i.e., that d{x)^/dm d — d(x)^/dm u and 
d(x} d /dm u — d(x)P/dmd, has been used to simplify the 
expression. 

Near the SU(3) symmetric point, the strange quark is 
considered as a 'heavy light quark', so that 



dm u 
9{xf d 
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mi 
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d{xf d 
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Rearranging, the CSV momentum fractions can be writ- 
ten as R7| 



1 m s ((a 



)/(x) 



Sd 



2 m 
1 ms 



q ( m 

((x)f 



-ml)/ 'X 2 
(x)f)/(x)\ 



2 m„ 



{m 2 K -ml)/Xl 



(12) 
(13) 



where m q — (2to/ + m s )/3 and X 2 = {2m 2 K + m 2 )/3. 



Similar expressions hold for the spin-dependent CSV mo- 
ments. This method allows an estimate of CSV at the 
SU(3) symmetric point. 

To evaluate the CSV terms at the physical point, we 
perform a chiral extrapolation of the lattice data for the 
quark moments 17]. As the isospin- averaged and -broken 
expressions for the Mellin moments as functions of pseu- 
doscalar or quark mass have the same unknown parame- 
ters, a fit to the available isospin-averaged lattice results 
allows the CSV terms to be evaluated from the isospin- 
broken expressions using Eqs. ([2| and ^ - a technique 
also used in [23]. These expressions can be evaluated 
at any pseudoscalar masses, in particular at the physical 
point. 



III. EXTRAPOLATION OF LATTICE DATA 

A. Fit to isospin-averaged lattice data 

In previous work, we described an isospin-averaged chi- 
ral perturbation theory fit to QCDSF/UKQCD Collab- 
oration lattice data for the first few Mellin moments of 



quark distributions [S]. Complete details of the fit for- 
mulae, fit parameters and method are given in Ref. |17j . 

In brief, chiral perturbation theory expansions, de- 
scribed in Ref. [17], were fit to QCDSF/UKQCD Col- 
laboration lattice data for the first spin-independent and 
zeroth and first spin-dependent moments. The fit func- 
tions include loop corrections and counterterms to lead- 
ing non-analytic order. In particular, the effect of chiral 
loops with both octet and decuplet baryon intermedi- 
ate states, as well as, for the spin-dependent moments, 
loops involving a transition between octet and decuplet 
baryons, are included. Tadpole diagrams and terms rep- 
resenting wavefunction renormalization are also consid- 
ered. 

The finite-range regularization scheme (FRR) is used 
to regularize the loop integrals. This technique, dis- 
cussed further in Rcfs. 18-20 , involves the introduction 
of a mass scale A through a regulator u{k) inserted into 
each integral expression. A is related to the scale be- 
yond which a formal expansion in powers of the Gold- 
stone boson masses breaks down (this scale is typically 
~ A/3 for a dipole). For this analysis, a dipole regula- 
tor u{k) = ^ A 2 A | _ A .2 ) and a regulator mass A = 1 GeV 
are chosen. This is based on a comparison of the nu- 
cleon's axial and induced pseudoscalar form factors [22] 
and the value of A deduced from a lattice analysis of nu- 
cleon magnetic moments |21j . All results are insensitive 
to this choice; choosing different regulator forms, for ex- 
ample monopole, Gaussian or sharp cutoff, and allowing 
A to vary by ±20% does not change the results of the 
analysis within the quoted uncertainties. 

Within the FRR framework, expressions for loops with 
octet intermediate states involve the integral: 



Tf 2^ 4 [°° k 4 U 2 (k) 

J(m z ) = - / dk — . 

3i (Vfc 2 +m 2 ) 3 



(14) 



with the finite-range regulator u{k) inserted. The nor- 
malization of J(m 2 ) has been defined so that the non- 
analytic part matches the common form of dimensionally 
regularized results, as J(m 2 ) — !• m 2 ln(m 2 //i 2 ). Loops 

with decuplet intermediate states may be written in an 
analogous way, in terms of 



Tf 2 ^ 4f J( fc 4 « 2 (fc) 

Jl(m '^3j( dk Wk i + mWkUm* + 6) 



(15) 



and 



Mm 2 , 8) - 



dk 



k A u 2 (k) 



3 Jo (Vk 2 + m 2 )(y/k 2 + m 2 + 5)< 



(16) 



which describe loops with one and two decuplet propaga- 
tors respectively. The tadpole contributions are written 
in terms of 

J T (m 2 ) = A 1^ dk k y {k \ (17) 



\Jk 2 + m 2 



3 



which has the same non-analytic structure as J, i.e., 



J T (m 2 



DR 



m ln(m //J ). To make comparison with DR 



expressions clear and to avoid absorbing loop terms into 
known parameters such as F and D, constant terms are 
subtracted by the integral replacement 



l(m^) ->• T{m$) = [2(m0) - 2(?7i0 = 0)] 



(18) 



where 2 represents any of the integrals in Eqs. ( 14 )-( 17 1 



The fit to the lattice results is performed by minimiz- 
ing the sum of x 2 f° r each set of moments. There are 
24 lattice data points available for each moment con- 
sidered [H [53] ■ The fit parameters, discussed in detail 
in Ref. [T7] and listed in Appendix [Aj are different for 
each moment. For the zeroth spin-dependent moment 
there are eight free parameters, while both the first spin- 



dependent and first spin-independent moment have nine 
fit parameters. 

[6j [7j taken from Ref. [T7] and located in 



show the quality of fit for each moment. 
= y/(2m 2 K + ml)/3 = 411 MeV is the sim- 
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Appendix 

Here X„ — v v*"^ 
ulation centre-of-mass of the pseudoscalar meson octet. 
Ratios of moments are shown and the X v normalization 
is taken for the figures so that they may be easily com- 
pared against published results [S] ■ The quality of fit is 
clearly acceptable in each case, with the x 2 /dof between 
0.6 and 0.9 for each moment. All x 2 values are less than 
one as the effect of correlations between the original lat- 
tice data could not be included. Figures [I] [2] and [3] show 
the fits to the data in a form suitable for the extraction 
of the CSV terms at the unphysical symmetric point by 
Eqs. (12) and (13 1. The full analysis, presented in the 
next section, includes an extrapolation to the physical 
point. 
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FIG. 1. Illustration of the fit to the zeroth spin-dependent moments - data from Ref. [8]. 
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FIG. 2. Illustration of the fit to the first spin-dependent moments - data from Ref. [8]. 



B. Evaluation of CSV terms tuting the best-fit parameters of the isospin-averaged fit 

As described in section [TTJ the CSV terms given in 
Eqs. ^ and (|3) may be evaluated by simply substi- 
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FIG. 3. Illustration of the fit to the first spin-independent moments - data from Ref. [8]. 



described in section III A into the full SU(3), isospin- terms. For example, SAu m may be expressed as a func- 



broken, perturbation theory expressions for the CSV tion of quark mass in the form |17j : 



where 



dm) 



(m) 
9 A 



16tt 2 / 2 



l - (-A,< ra) + An { ™ } + Ant } + A"!™') B(m u - m d ), 



1 



6^3 
l r 

+ 24 



D 2 - 2DF - 3F 2 ) sin(2e) (5Aa (m) + 2A(3 ( - m '> + 6Aa^ J(ml o ) - J(m 2 ) 
D 2 (9Aa (m) + 2A^ m > + 8A<j( m >) + 2DF (l9Aa (m) + lOAfi^ + 24Act< 



3F 2 ( 5Aa (m) + 2A(3^ + 8Aa (r ' 



J(m 2 K o) - J(m 2 K ±) 



+ - (5Aa^+2A^) [j T (m 2 K0 ) - J T {m 2 K± ) 



d { r ] = 



i 

72 
1 

108 

4 
9%/3 



(5Aa^ + 2A/3^ + 6Aa^) C 2 [j 2 (m|o, <5) - J 2 {m 2 K± ,5) 



5 A 7 W _ A 7 '( m )j C 2 [j 2 (m 2 K „,S) - J 2 (m 2 K± ,S) 
(D - 3F) sin(2e)A W (™) [j^m^o, S) + Ji{m 2 5) 



+ -(D- 3F)AljW J x (m 2 K0 ,$)-Jx (m 2 K± , S) 
y 



(19) 

(20) 
(21) 



(22) 



(23) 



and expressions for the (subtracted) integrals J are 
given in the previous section. Clearly, entirely analo- 
gous expressions may be written for 5Ad m and the spin- 
independent CSV terms. These, taken from Ref. [17], 
are given in appendix |Aj We remind the reader that, to 
the same order in the broken SU(3) symmetry, analo- 
gous expressions for each quark flavour combination in 



each octet baryon are expressed in terms of different lin- 
ear combinations of the same coefficients. The general 
case is given in Ref. [T7| . 
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In the above expression, meson masses take the form 

(24) 



values at the physical point at 4 GeV 2 



m 2 K± 
m 2 K o 



-B(m u 
--B(m u 
2B 
~ ~3~ 
-B{m a 
--B(m s 
B 
~3 

2B 

+ —(2m s - (m u + m d ) ) 



+ m d ) 
+ m d ) 

(2m s - 

+ m u ) 
+ m d ) 



-(4m s 



(m u + m d )) 



m d ) 



■ 2 

sin e 



cos2e 



sin 2 e 

3 cosze 
and the ir° — rj mixing angle e is given by 

(m d - m u ) 



tan2e 



2m s - (m d + m u ) 



(25) 

(26) 
(27) 

(28) 
(29) 



The parameters An^ , AaW, A/3^ and Aa^ are 
determined, for m = {0, 1}, from the isospin-averaged 
fits. All that remains to be specified for an evaluation of 
SAu m from the expression above are values for Bm q . 

To evaluate the CSV terms at the physical point we 
take as input the estimate for the physical up-down quark 
mass ratio from Ref. [25] 



R 



m d 



0.553 ±0.043, 



(30) 



determined by a fit to meson decay rates. We note that 
this value is compatible with more recent estimates of the 
ratio from 2 + 1 and 3 flavor QCD and QED [26][27]. and 
lies within uncertainties of the FLAG lattice averaging 
group estimate R = 0.47(4) [28 . The Gell-Mann-Oakes 
Renner relation suggests the definition 



B(m d - m u ) 



which allows one to define 



iii-^m 2 (31) 



Bm u 


= ml 12 

^(phys)' 


- u, 


(32) 


Bm d 


= ml 12 

T(phy a ) ' 


+ w, 


(33) 


Bm s 


= ml- — 

«( P hys) 


ml 12. 

lT(phy a ) ' 


(34) 



Here, m^^, = 137.3 MeV and m K{phys) = 497.5 MeV 
are taken to be the physical isospin-averaged meson 
masses [29] . 

As the available QCDSF/UKQCD Collaboration lat- 
tice results are presented only in terms of ratios of mo- 
ments, there is an unknown constant scaling factor Z as- 
sociated with all data points. This Z is distinct for each 
moment (zeroth and first SD and first SI) under consider- 
ation. These constants are determined by matching the 
extrapolations for the isovector moments to experimental 



9A = (l) P Au - Ad = pt 1-2695(29), 
^u-A d = t °- 19 °( 8 )' 
^^ cx = pt °- 157 ( 9 )- 



(35) 
(36) 
(37) 



The uncertainty of these experimental numbers is 
propagated into the final results. The full error analy- 
sis also takes account of correlated uncertainties between 
all of the fit parameters in the original fits [17], as well 
as allowing for the stated variation of R. The regulator 
mass A = 1 GeV is allowed to vary by ±20%, which is 
again propagated into the final uncertainty. Changing 
the regulator u(k) within the FRR scheme leads to small 
variations of order 1%. 

The results of this analysis are summarized and 
compared with previous work in Table [T] While the 
light quark ratio R was used as input in this cal- 
culation, the determination of the CSV terms via a 
linear flavour expansion [5] used the quark mass ra- 
tio 3(m d — m u )/(m d + m u + m s ) = 0.066(7) 25J. The 
choice of R made here sets this ratio to the same value. 
The other inputs used in both calculations, namely the 
experimental isovector moments at the physical point, 
take the same values in both calculations. Thus, the lin- 
ear and chiral results in Table [I] are directly comparable. 

In particular, evaluating the chiral perturbation the- 
ory expressions for the CSV terms at the point where 
(m d +m u ) = 2m s and both {m d —m u ) and (m u +m d +m s ) 
take their physical values, labelled 'SU(3)-sym' in Ta- 
ble [TJ gives results which may be directly compared with 
the linear flavor expansion calculation. As might be an- 
ticipated from an inspection of Figs. [TJ|3] which show fits 
qualitatively consistent with straight lines, chiral loop 
corrections to the CSV terms at this point are small and 
within uncertainties. 

Comparison with results evaluated at the physical 
pseudoscalar masses gives an indication of the chiral loop 
corrections in moving away from the SU(3) point. Again, 
these corrections are small in the spin-dependent case, 
while being more significant in the spin-independent case. 
It is noted that, in contrast to the results of the linear 
flavour expansion, the chiral perturbation theory results 
are based on fits for the quark distribution moments to all 
lattice data simultaneously (for each moment), and thus 
include the proper correlations between quark moments 
in each of the baryons. As a consequence, even with more 
fit parameters, the uncertainties are comparable to the 
simple linear fits. 

The origin of the chiral loop contributions to the CSV 
terms can be seen clearly from the form of Eq. ( 19 ) (and 
the analogous Eqs. ( Al ), ( |A6[ ) and ( A7) in Appendix [A|. 
One contribution to the [u — d) moments is illustrated 
diagrammatically in Fig. [4j The kaon loop diagrams 
shown, and the analogous diagrams for the (d — u) mo- 
ments, give contributions to the CSV terms proportional 



to 



J(m 2 K0 ) 



J(m 2 K ±) 



which is non-vanishing when 



() 



Moment 



5Ad°+ 

SAu 1 - 

SAd 1 - 

5u 1+ 

5d 1+ 



Linear: SU(3)-sym 



Chiral: SU(3)-sym 



Chiral: physical 



-0.0057(14) 
-0.0018(6) 
-0.0010(3) 
-0.0004(1) 
-0.0012(3) 
0.0010(2) 



-0.0063(13) 
-0.0019(6) 
-0.0007(2) 
-0.0003(1) 
-0.0013(3) 
0.0012(2) 



-0.0061(13) 
-0.0018(6) 
-0.0007(2) 
-0.0002(1) 
-0.0023(7) 
0.0017(4) 



TABLE I. Comparison of results. The column labelled 'Linear' gives the results which were published with the lattice simulation 
results [B], calculated using a linear flavor expansion about the SU(3) symmetric point. These have been corrected from the 
values quoted in [8j, as explained in the footnote preceding Eq. (12 I. 'Chiral' gives the results of this work, i.e., including 
chiral physics, both at the comparable 'SU(3) symmetric' point (with (md + rn u ) = 2m 3 but the physical (m<j — m u )), labelled 
'SU(3)-sym', and at physical pseudoscalar masses. 



m K o m K ± . The corresponding wavefunction renor- 
malization terms, as well as tadpole and decuplet kaon- 
loop diagrams, also contribute to the CSV terms propor- 



tional to 



J{m 2 K0 ) 



J(m 2 K± ) 



. In the spin-independent 



case, this kaon mass difference effect yields the only chi- 
ral loop corrections to the CSV terms. For the spin- 
dependent moments, however, additional terms propor- 



tional to 



also contribute. Cancella- 



tion of octet loop terms with wavefunction renormaliza- 
tion contributions ensures that these terms vanish in the 
SI case. 

The chiral loops also account for the corrections in 
moving from the l SU(3) point' to the physical point. For 
example, as one moves along the line of constant sin- 
glet quark mass {{m u + nid + m 8 ) = constant) from the 
SU(3) symmetric point to the physical point, the dif- 



ference 



approximately 30%. 



J{m 2 K0 ) — J(m 2 K ±) decreases in magnitude by 



A' + 



K" 



, u 

-m x- 



A,E° A,E° 



d 

-x- 



A.S° A,S° 



A 



u 

-x- 



d 

-x- 



FIG. 4. Illustration of some of the octet loop terms contributing to SAu m = {x m ) p ^ u - 
contributions are non- vanishing when the loop pseudoscalar masses are different, i.e.. 



when m 



x m )<X, or 8u m = (x m )l- 



2 _L ™2 
K° 



{x m )%. These 



m K± . 



r 



IV. CONCLUSION 

We have used a chiral perturbation theory analysis to 
extrapolate QCDSF/UKQCD Collaboration lattice data 
for the first several Mellin moments of quark distribution 
functions to the physical quark masses. This technique 
allows the charge symmetry violating (CSV) parton dis- 
tributions to be evaluated at the physical point. 

The conclusion of this study is quite clear. The chiral 



corrections to the spin-dependent CSV moments are very 
small. In particular, an analysis of the same lattice data 
using a linear flavor expansion about the SU(3) symmet- 
ric point gave compatible results [5]. A detailed analysis 
shows that both the chiral corrections at the SU(3) sym- 
metric point, as well as the extrapolation from this point 
to the physical quark or pseudoscalar masses, are small 
effects. 

At the physical point, this analysis gives the spin- 
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dependent CSV terms to be 5Au°+ = -0.0061(13), 
8Ad 0+ = -0.0018(6), SAu 1 - = -0.0007(2), and 
6 Ad 1 - = -0.0002(1). As a result, one would expect CSV 
corrections to the Bjorken sum rule [33J 133] to appear 
at the half-percent level. Measuring these corrections 
would require significant improvement over the current 
best determination of the sum rule to 8% precision at 
Q 2 = 3 GeV 2 from a recent COMPASS Collaboration 
experiment |35) . 

For the spin-independent moments, the chiral correc- 
tions are more significant. This analysis gives Su 1+ = 
-0.0023(7) and 5d 1+ = 0.0017(4), in good agreement 
with previous phenomenological estimates of CSV both 
within the MIT bag model |7] and using the MRST 
analysis [3] . These results support the conclusion [3J [7] 
that partonic CSV effects may reduce the 3a discrepancy 



with the standard model reported by the NuTeV Collab- 
oration [36] by up to 30%. 
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Appendix A: Extrapolation formulae 



This section gives formulae for the spin-dependent and spin-independent charge symmetry violating quark distri- 
butions as functions of quark and meson mass. All integrals are defined in the body of the report. 

The fit parameters which appear in the following expressions are discussed in detail in Ref. [T7]- For the zeroth 
spin-dependent moment there are eight free parameters; six linearly independent linear coefficients An-^.g, the 
baryon-baryon-meson coupling constant D, and an operator insertion parameter Act' '. Both the first spin-dependent 
and first spin-independent moment have nine fit parameters; six linear coefficients An^ 1 _ 6 (ri|2 1 _ 6 ), and three 
operator insertion parameters Aa^\ A/3^\ Aa^ (a^ 2 \ fl( 2 ',<j( 2 >) in the spin-dependent (-independent) case. 

Baryon-baron-meson couplings F and D are set to their physical values by D — > D p h ys = |gA phya for each of the 
first-moment fits. For all three fits, SU(6) symmetry is used to set F — |D and C — > C p h ys = — §3A phys is fixed. 
Decuplet (7) and transition (w) insertion parameters are also fixed for each fit, either by using SU(6) symmetry to 
relate them to other fit parameters, or, in the case of 7 for the first spin-independent moment, by using an experimental 
result, as detailed in Ref. [P7] . 



1. Spin-dependent CSV terms 



This section gives an explicit expression for the spin-dependent CSV distribution SAd m as a function of quark and 
meson mass. The corresponding expression for 8Au m is given in the body of the report. 



5 Ad™ = (x m )l d - (x™)l u = + _L_ (4 m) + + (Al) 



_(m) 

1 v ' 



i(-A4 m) +A4 m) )B(m u -m d ) 



r(m) 1 



6^ 

1 

+ 24 



(D 2 - 2DF - 3F 2 ) sin(2e) (Aa< m ) + 4A/3( m ) + 6Acr(" 1 )) [j(m 2 J - J(m^) 

L> 2 (A«< m ) - 4A/3< m ) - 8Acr< m )) + 6DF (Aa'™' + 4A/?( m > + 8A<r^ 
+F 2 (5Aa (m) + 20A/3 (m) + 24Aa (m) )] [j(m^ ) - J(m 2 K± ) 
- 1 (A a ( m ) +4A/^ m )) [j T (m^,) - Jr(m 2 K± )" 



di m) = - — (Aa^ +4A/3< m ) +6Aa( m )) C 2 [j 2 (mW) - J 2 (m 2 K± ,S) 
+ (5A 7 (™) - A 7 '(™)) C 2 [J 2 (m| ,5) - J 2 (m^ ± ,<f) 



4 

9^ 

4 



3A m) = + TT/f (£> - 3F) sin(2e)Ao/ m ) Ji(m 2 o, <5) + J^m 2 , 5) 



9 



J\{m 2 K0 , 8) - Ji(m 2 K ± , <5) 



(A2) 
(A3) 



(A4) 



(A5) 



2. Spin-independent CSV terms 



This section gives explicit expressions for the spin-independent CSV distributions as functions of quark and meson 
mass. 



6vT 
Sd n 



-{x m )l - (x m ) n d = « (m) + ^j- 2 (b {m) + d (m) ) 

__ {xm) p „ {x m } n = n (m) + _J_ + ^ 



(A6) 
(A7) 



1 \ 1 / 

» W = 2 


-»i 






24 






l / 


^24V 



J{m 2 K a) - J(m 2 K ±) 



5a 



(™) + 2/?( m >) [j T (m 2 K o) - J T (m 2 K ±) 



36 



^3 7 (m) _ y(m)) C 2 [j^^) _ J 2 ( TO ^ ± , £) 



b v " =— 



24 



1 

24 



(a (m) + 4/?(" 1 )) [j T (m^o) - J T (m 2 K± ) 



- ^ (a< m > +4/3 (m) +6<r< m) )c 2 [j 2 (m^o,<y)- J 2 (m^ ± ,<5) 
+ ^ (3 7 (m) - 7 M ) C 2 [J 2 (m|a,5) - J 2 (m| ±I <5) 



9 

(A8) 
(A9) 

(A10) 



%) =\ (-"3 m) + S(m„ - m d ) (All) 
" ' f -D 2 (ja^ + 4/3<" 1 )) + 6DF (a™ + 4p (m) ^) - 3F 2 (a^ + 4/3^ J(m 2 K „) - J(m 2 K± ) (A12) 



(A13) 



Appendix B: Figures 

This section shows the fit s to Q CDSF/UKQCD lattice 
results discussed in section |III A| The figures are taken 



from Ref . |T7j , and are included here to give an indication 
of the quality of fit. 
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